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^ . Abstract 

| The interaction potential between the alpha particle and the de- 

O ' formed parent nucleus was used for description of the decay of super- 
heavy nuclei. It consists of centrifugal, nuclear and Coulomb parts 

| suitably modified for deformed nuclei. The significant effect of various 

^> ■ shapes of barriers obtained for deformed parent nuclei on calculated al- 

, pha half-life times were shown. The finally calculated half-life times due 

^T* ! to the spontaneous alpha decay for superheavy elements were compared 

q ■ with the results obtained from other models and the experimental data. 

3 : 

C ■ 

1 Introduction 

X' 

Theoretical calculations suggesting existence of deformed superheavy nuclei 
from the areas of doubly magic 270 Hs were made by Sobiczewski et. al.[l, 2]. 
Recently there have been calculated some collective properties of nuclei such as: 
deformation energies with the equilibrium shapes, energies of the first excited 
state E2+ and the branching ratios between a-decay to this state 2+, to that to 
the ground state 0+. Possible puzzling comparison of the obtained theoretical 
predictions with the experimental data promotes and motivates further studies 
of this problem. Therefore more extensive description of an a particle emission 
including additional degrees of freedom connected with deformation in ground 
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state of emitting nucleus is necessary. Then different possition of a particle 
towards this nucleus is possible. 

It is commonly known that the action integral giving probability of a par- 
ticle tunneling depends strongly, among others, on the potential barrier shapes 
and even insignificant change of this barrier can change the obtained half-lives 
by the magnitude of the a few orders. Taking deformation of parent nucleus 
into account in our model can give more realistic results. 

The aim of the paper is calculation of half-lives due to spontaneous a de- 
cay for the deformed superheavy nuclei which is a predominant process for 
them. The barriers were calculated using a suitable model in which the total 
interaction potential (Vtot) consists of the nuclear (Vn), Coulomb (Vc) and 
centrifugal barrier (Vj,) modified for the case of existing the deformation of 
parent nucleus. Chapter II includes brief description of all terms of the poten- 
tial. Chapter III presents main formulas relating to the a decay used in this 
paper. Chapter IV includes the discussion of the results and the influence on 
half-lives of the various factors. The last chapter includes summing up and 
conlusions. 



2 Interaction potential 

2.1 Shape parameterization in the exit channel 

The shape of nucleus surface area -Rp( a )(0p( a ), 4>p(a)) is parameterized in the 
standard way in the spherical harmonic base which can be written in the vector 
notation as follows: 

!Rp( a ){6p( a ), <t>p{a)) sin(0 P(a )) COS(0P( Q )) 
Rp(a)(0p(a), 4>P(a)) Sm(0p( a )) sin(0 P ( a) ) (1) 
Rp(a)(0p(a), (f>P(a)) COs(0p( Q )) 

The index P(a) refers to parent nucleus or alpha particle adequately. In the 
reference system connected with the parent nucleus the centre of the emitted 
particle is parameterized in the spherical system (R, 0, $): 

f J Rsin(0)cos($) 
R= i #sin(0)sin($) (2) 
( J Rcos(0) 

Figure 1 shows geometry of this system where the shortest distance between 
the nucleus surface areas of a particle and the parent nucleus is denoted by 
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s, the distance between the mass centers by R and the distance between the 
infinitely small volume elements in the nucleus of emitter and emitted particle 
by r Pa . 

Particle a is treated as a sphere of a radius R a = 1.84/m. There were 
considered only axially symmetric deformed parent nuclei which results in the 
simplification - 4>p( a ) — and $ = 0. 



2.2 Coulomb potential 

The electrostatic interaction potential (Vc) between the parent nucleus with 
the proton (Zp), mass numbers (Ap) and the emitted a particle is given by 
integral equation: 

Vc = 2e[ P -^ldT P , (3) 
J up r Pa 

where e is the charge unit and pp(fp) is a charge distribution with uniform 
density It is pp(r P ) = inside the parent nucleus and p P (r P ) = outside 

CP ^ 

the nucleus. The equivalent sharp radius is taken as Rep = 1. 15 Ap. Integra- 
tion over the parent nucleus volume (up) in the coordinates defined by equation 
(1) is carried out. In the outside area, when the distance between the mass 
centers is larger than the distances corresponding to the contact configuration 
of the parent nucleus and the alpha particle R > r P + r a , expansion of distance 
reciprocity r Pa is of absolutely convergent series form: 

After substitution of the expansion into the formula (3) and making simple 
transformations we obtain suitable multipole expansion for the Coulomb in- 
teractions of deformed nuclei: 

Here C£ 'o are the Wigner coefficients. The required multipole moments Q® 
are calculated using the following definition [6]: 

Ql = I Pp{?pV?Y? P {0 P )dTp . (6) 

J Up 

here Y l ° p (6 P ) are corresponding sherical harmonics. Let us note that the first 
term in the presented multipole expansion corresponds to the case of interac- 
tion when both fragments are spherical. 
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2.3 Nuclear potential 

The nuclear part (Vat) was calculated using proximity aproximation [7]. This 
hypothesis is based on the assumption, that the nuclear interaction takes place 
between two infinitesimal volume elements located closest for a given configu- 
ration, and additionally that the distance between them is of the order of a few 
fm. The surfaces could be treated then as flat parrallel slabs and the proximity 
energy between two semi-infinite nuclear matters is of a simpler form: 

V N = K 1 *P(s) , (7) 

with the surface tension coefficient 7 = 0.9517(1 — 1.7826I 2 )MeV/ fm 2 , where 
I — (N — Z)/A. It allows for division of nucleus- nucleus interactions into two 
independent factors: K connected with the shape of emitting nucleus surface 
area and the factor expressed by the universal function ip(s), independent of 
the system geometry The form of universal function given in [8] is used in 
this paper. 

Expansion of this approach for the deformed case consists in suitable re- 
defining of the geometrical factor. In the case of deformed nuclei it is the 
geometrical mean of main reduced curvature radii k\ and k 2 : 

K DEF = 2n^k^k~ 2 , (8) 

calculated on the nucleus surface area in the sites for which with the data 
R i G the distance between surfaces is minimal s = min\R + f a — fp\. To 
calculate the required curvature radii, we calculate fundamental forms of the 
first and the second orders of the nucleus surface in these points. Then we 
are able to calculate the Gauss curvature and the mean curvature and finally 
we obtain the first Rim (a) an d second R2M( a ) main radius curvatures [9]. The 
main reduced curvature radius of the first and second orders k\ and k 2 can be 
expressed by the relations: 

kl = r — Tr~ ] h = r — Tr~ ' (9) 

Moreover, within small deformations the factor Kdef passes into the factor 
K for the spherical case K DEF (f3i p — 0)— > K. 

3 Theory of alpha decay 

The emission of alpha particle from the axially symmetric deformed superheavy 
nuclei is considered. The half-life is given as a function of the position a 
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particle: 



71/2(0) 



ln2 



(10) 



A(0) ' 



where A is the decay constant. The a-decay model for deformed superheavy 
elements presented here is based on the assumption that forming of a cluster 
in the inside of parent nucleus and passing over potential barrier are inde- 
pendent events. It is also assumed that the emission of a particle does not 
affect deformation of the nucleus from which it was emitted. Within these 
approximations the a-decay tunnelling probability constant can be writen as: 



where S is the quantum mechanical probability forming of a particle in the 
parent nucleus (preformation probability) . It is additionally assumed that the 
deformation and orientation of parent nucleus do not have any influence on the 
preformation probability. By the Ag(O) we denoted the tunneling probability 
(Gamow decay constant), which is a function of the energy of the a particle 
and angle showing the direction of a emission. All informnation about the 
many body structure is contained in the S and because of it we sometimes 
call this constant the spectroscopic factor. This factor can be calculated as the 
mean value of the suitable projection operator P between parent nucleus wave 
function $p lies in the open channel space $£> <S> spread by the daughter 
nucleus wave function $£> and a nucleus wave function $ a . 



This problem was invesigated by Blendowske, Fliessbach and Walliser [12] for 
different exotic spontaneous emissions of clusters and for a particle needed 
here. Exact numerical evaluation and fitting procedure lead to the semiempir- 
ical formula for spectroscopic factor. In the a particle case this factor is equal 
to S a = 6.3 • 10~ 3 [12]. The favoured decays are considered here only. The 
barrier penetration probability \g as a funcion of the angle can be expressed: 



A(0) = SA G (0) , 



(11) 



S = ($p|P|$ P ) . 



(12) 



A G (0) 



<4®) 7(G) 



(13) 



where J(0) is the action integral: 




(V TOT (R,Q)-Q a ) dr, 



(14) 
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The reduced mass of a particle is denoted by /i, r en (0) and r ex (Q) are 
the classical turning points determined for a given position of a particle by 
the relation Vtot{R, 9) — Qa- The energy released when the nucleus (Z,N) 
emits an a particle to the state L+ is obtained from the mass or total binding 
energies: 

Q a (Z, N) = E B (Z -2,N-2)- E B (2, 2) - E B (Z, N) - E L+ (Z, N) , (15) 

where El + is the rotational energy of the L+ state of the nucleus (Z,N). The 
total potential (Vtot) consists of Coulomb (Vc), nuclear (Vjv) and centrifugal 
(Vl) energies: 

V TOT (R, 9) = V C (R, 9) + V N (R, 9) + V L (R, 9) , (16) 
The centrifugal term of energy can be calculated: 

v(Re) h 2 L(L+l) 

Vl ^°>~ 2jRW 2 ' ( } 

The knocking frequency is a slowly varying function of E. 



u < e > = *VOTr (18) 

In our investigations we accepted the a particle kinetic energy E = 112.1 MeV 
and a particle binding energies E B (2,2) = —28.30 MeV according to [10, 11]. 
The obtained half-lives were averaged with the surface function in order to 
compare them to the experimental data. This surface function in the axially 
symmetric geometry case is given by: 



2 

s = 2,/ orfS( e) = 2. j ( de^ejsm^i + ^^J , (in) 



The mean half-life can be writen in the following way: 

£ T 1/2 (&) dS (Q) 



4 Results 

As follows from the formulae presented in the previous chapter, the precise 
knowledge of the barriers shapes depending on the angle 9 at which the par- 
ticle is emitted from the nucleus is necessary for the correct estimation of the 
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tunnelling probability (formulae 13,14). Figure 2 shown the obtained total bar- 
riers Vtot{R-, O) in the function of the distance between the mass centres of a 
particle and the emitter for the deformed of the superheavy element 264 Hs^q. 
The equilibrium deformations {/?2, Ai, (3q} became taken from [14]. The barri- 
ers for different orientations 6 of the emitted a particle in the relation to the 
parent nucleus are shown. Note that position of the barrier maximum and its 
value change with the angle of a cluster emission as shown in the upper inside 
panel in figure 2. The lowest barrier height in the given example appears at 
the angle 6 ~ 30° but the highest one is obtained for the equatorial position 
6 90°. The difference in the barrier maxima is AVmax = 1.7 MeV for this 
nucleus. It has significant influence on the half-lives as the difference in the 
barrier height being only AVmax — 0.5 MeV gives approximately the half-life 
differing the order of magnitude. 

Another parameter to which emission probability is very sensitive is the 
entrance energy of a particle in the barrier Vent — Vst + Qa where Q a is the 
a-decay energy calculated from the formula (15) and is equal Q a = 10.57MeV 
for the studied isotope 264 Hs\^. Vst is the energy of potential corresponding 
to the configuration in the contact point between a particle and the parent 
nucleus {sticking configuration). It is functional dependence on the angle 9 
is shown in the bottom panel in Figure 2. Taking into account the presented 
dependences leads to angular anisotropy and the largest difference is AVent = 
2.5 MeV. As previously it can lead to essential differences in the values of 
action integral, the obtained tunnelling probabilities and further in calculated 
half-lives. 

In order to show the effect of the deformation parameter /3 6 on the obtained 
half-lives the calculation for both cases /3 6 > and /3 6 < (|/3 6 | = 0.034) was 
made leaving the deformation parameters {A2, Ai} unchanged. Figure 3 show 
the course of variability of tunneling probability function 1(0) in both cases. 
Analyzing the diagrams it can be seen that position of the curve maximum 
at (5q < denoted with the solid line correspond to the position of minimum 
for the curve at f3 e > denoted with the broken line. Position of the curve 
maximum for positive (3q sign is shifted by about 15° towards smaller angles 
compared to the curve minimum with negative sign of f3§. As follows from 
the above analysis half-lives are very sensitive to even smallest changes of 
deformation parameter /3 6 . 

On figure 4 calculated half-lives in the function of angle 6 denoted by 
the thick solid line are shown. The value (7i/ 2 ) obtained after averaging in 
relation to the surface according to (20) is also denoted (thin solid line). The 
arrow indicates the experimental values [16]. Half-lives differ from each other 
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depending on the place of a particle emission even by 5 orders of magnitude and 
the obtained mean value of half-life for this element is close to the experimental 
value. 

The dotted line correspond to the calculus without taking centrifugal po- 
tential into consideration (emission from ground to the ground state is con- 
sidered). The calculated energy of the 2+ state for this element is very small 
beeing E 2+ = 46.6 keV see eg. [3]. One can also find that the contribution of 
the centrifugal barrier is small and it can be neglected. 

Systematic calculus of half-life for superheavy elements Z=108,110,112,114 
is presented in figure 5 where decimal logarithms of the average half-life (dark 
filled squares) in the function of the mass number A are shown. The obtained 
results are compared with literature calculations based on the phenomenologi- 
cal formulae given by Viola and Seaborg (V-S) [13, 15] denoted by crosses, on 
the Generalized Liquid Drop Model (GLDM) [17] marked with blank squares 
and with available experimental values (EXP) [16] given by full triangles. The 
calculations were also made not taking deformation into account denoted with 
open circles. As can been seen not taking deformation into account leads to 
the lowest half-lives while including deformation degrees of freedom increases 
significantly the values of half-life. The best agreement with the experimen- 
tal data is observed for the isotope with Z=110. For the element Z=108 the 
obtained results are the closest to the predictions obtained using the Viola- 
Seaborg formulae [13, 15] and for the element Z=112 the largest agrement 
was obtained with calculations based on the GLDM model. The comparison 
of our results and those obtained from the Viola and Seaborg model shows 
that the latter are systematicaly higher. The results for the isotope Z=114 
obtained by us lie below the values obtained using other models. Thus for the 
lately synthesized element 289 114 the phenomenological Viola-Seaborg formu- 
lae gives the value of half-life time logarithm log(T\j 2 ) = 4.86 s while that cal- 
culated in this paper is logiTi^) = —1.30 s and the experimental data indicate 
log(Ti/ 2 ) = 1.32 s. For the neighbouring isotope 288 114 the experimentally pre- 
dicted half-life is log(Ti/ 2 ) = 0.26 s see [18] ( and citations in it) but we obtain 
the value almost three orders of magnitude smaller log(Ti/ 2 ) = —2.86 s and es- 
timation based on the Viola-Seaborg formulae gives the value log(T\/ 2 ) = 2.80 s 
almost three times as large. This can be explained by small equilibrium defor- 
mation of superheavy nuclei under discussion. The results obtained for them 
are similar to those obtained for the spherical case and as mentioned earlier 
they are usually lower than the actual ones. Despite some discrepancies the 
values of the half-time obtained by us not based on any phenomenological 
formula are similar to the experimental values. 
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5 Conclusions 



The following conclusions can be drawn from our investigation devoted to a- 
decay of deformed superheavy elements: 

1. Taking deformation of nuclei into account affects significantly the poten- 
tial barrier shapes: position of maximum and entrance energies and in 
consequence, the half-life time. 

2. Even small changes of deformation (3q modify greatly the values of action 
integral. 

3. The half-life time for superheavy elements of Z=108-i- 114 without appli- 
cation of any phenomenological formula but taking into account defor- 
mation of nuclei in the ground state reproduce the experimental values 
well. 

4. The centrifugal term affects the obtained results insignificantly. 
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Figure 2: The total potential (Vtot) energy barriers obtained within the pre- 
sented model (solid lines) dependent on the distance R between the centers of 
the a particle and the daughter nucleus for different orientations. The upper 
panel correspond to the maximum of the barrier (Vmax) depending on O, The 
bottom panel gives the entrance energy (Vent) during the emission process. 
All presented results are predicted for the exemplary deformed superheavy 
264 if Si 56 nucleus. 
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Figure 3: The influence of emitted nucleus curvature on the obtained results 
is presented. The penetration probability as a function of the position of a 
particle for curvature with @ e = —0.034 (solid lines) and with /3 6 = 0.034 
(dashed lines) for the deformed 264 ifsi 56 is shown. 
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Figure 4: The logarithm of the a-decay half-life as an orientation function 
(thick lines) is shown. Solid lines indicate that the centrifugal barier is included 
in calculation but the dotted ones that it is not. The mean half-life (7i/ 2 ) is 
marked by the thin lines. The solid line is related to the calculation with the 
centrifugal term and the doted ones to that without this term. Comparison 
with the experimental data (EXP) is indicated by the arrow [16]. 
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Figure 5: Logarithm of the a-decay mean half-lives of the isotopes with atomic 
number Z = 108, 110, 112, 114 as a function of the mass number A. Full squares 
indicate the calculated and averaged half-lives (T1/2) • Open circles are related 
to our calculation without deformation (SPHERICAL), cross symbols denote 
the calculation with the semi phenomenological Viola-Seaborg (V-S) formula 
[13] with a new set of parameters given in [15]. The open squares were cal- 
culated using the GLDM model [17]. The triangles indicate the experimental 
data [16] (EXP). 
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